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Abstract 

We study parametric resonant decay of inflaton field with emphasis on 
its physical meaning. We show that the parametric resonance is indeed an 
induced process, which means that the more numbers of produced particles, 
the more inflaton field decays. We also consider the dissipative effects of 
produced particles and find that the dissipation reduces the resonant decay 
rate of inflaton field. 



1 Introduction 



Inflation |], [| is the most successful theory to solve many problems of the standard 
"hot big bang" Universe. In the inflationary Universe scenario the universe is ex- 
panded exponentially by the vacuum energy of some scalar field <fi (inflaton). For 
this scenario to work, the radiation-dominated Universe must be recovered after in- 
flationary stage. Therefore, some process is necessary to transfer the vacuum energy 
to relativistic particles. This process is called reheating. 

The old version of reheating theory was first considered in || for the new 
inflationary theory ||, and it can be applied to the chaotic inflation ||. After 
inflation, inflaton field oscillates near the bottom of its effective potential. Inflaton 
field decays into other (lighter) particles due to its coupling to others. The reheating 



temperature can be estimated as Trh — 10 _1 yr tot M p J7|, based on the single- 
particle decay. However, recent investigation revealed that the drastic decay of 
inflaton field occurs in the first stage of reheating || ||, [H], [11], [12], This stage 
is called preheating ||. At first, inflaton field (j) explosively decays into some bosons 
X, whose spectrum is far from equilibrium. After x-particles collide with each others 
or decay into lighter particles, thermal equilibrium can be achieved and reheating is 
complete. 

The equation for the particles produced can be described by Mathieu type equa- 
tion, which has instability solution in some regions of parameters (instability bands). 
If the relevant parameters stay in the instability bands long enough, the solution 
will explosively grow, which means that the number of particles created becomes 
exponentially large so that the parametric resonant decay takes place very efn- 



ciently |, |, O 11, 



From the different viewpoint, the physical meaning of parametric resonant decay 
could be considered as an "induced" decay, which is similar to an induced emission 
of photons in the two energy level system of atoms, in the sense that the presence of 
produced particles stimulates inflaton field to decay into those particles. Thus this 
phenomenon is peculiar to those particles that obey Bose-Einstein statistics, bosons, 
the number of which in each mode k can increase exponentially. This is the reason 
that parametric resonant decay into fermions cannot occur due to Pauli's exclusion 
principle 0, §, |, [TDJ . 



If the parametric resonance is an induced process, the parametric resonant decay 



1 



of inflaton into bosons are likely to be suppressed by the processes of scattering or 
decay of the produced particles [14], 3. The reason is that dissipation (the particle 
scattering, decay and annihilation) will tend to decrease the mean occupation num- 
bers Nk in each particular mode so that the induced decay will be reduced, and, as 
a result, these numbers cannot grow so large during the resonance period. 

Therefore, in the present paper we study the parametric resonant decay in the 
narrow resonance region [] with emphasis on its physical meaning and consider the 
effect of dissipation. In section 2, we will make clear that the parametric resonant 
decay is an induced effect, and we will see the restriction to that decay due to effects 
of dissipation in section 3. Section 4 contains our conclusion and discussion. 

2 Parametric Resonant Decay as an Induced Ef- 
fect 

In this paper we assume the interaction between inflaton field and x-particles as 
C-int = —vfix 2 an d the effective potential of inflaton field as V(<f>) = |m^0 2 . The 
similar results can be obtained for other interactions and also for potential. The 
mode equation for x held becomes 

Xk + K + 2or 0o sin(m t)]xfc = 0, (1) 

where oj\ = k 2 + m 2 , with m x being the mass of the x-particle, 0o is the amplitude 
of the inflaton field, and the expansion of the Universe is neglected. 
We can expand x field as 

X = j^j~ 2 J d 3 k[a(k)x4t) e ife + a\k) X+ (t) e"^], (2) 

where a?(k) and a(k) are creation and annihilation operators, respectively, with 
commutation relation: [a(k), a'(k')] = 5(k — k'). x± are the positive and negative 
frequency solutions of Eq. ([[]), which are equal to the positive and negative frequency 
asymptotic solutions (= exp(±k t) / y/2k~o) if x was free field. 

Using Eq. (fj) the energy density of the x-particle can be calculated as 

<o|p x |o> = (0 1 \x 2 + \^lx 2 1 0) 



f2vr) 3 



2' v 2 



d 3 k 
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(3) 



^^We study onlyin the narrow resonance region since we expect that the suppresion will be 
unimportant in the broad resonance region. 
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where we take the vacuum expectation value. Thus the energy density for the x~ 
particle with momentum k is given by 

pi = \\x-? + \^l\xA 2 - (4) 

Then the time evolution of p k x is described as 

^-rftU 1 ^ 1 J~2 dt lX ~ l ' (5) 

where the Eq. (D) is used in the second step. 

Hereafter, we assume that o~</> <C m| which corresponds to so-called narrow 
resonance region. This assumption allows us to take perturbative analysis. f] First, 
we will show that our analysis is general so that it does not depent on the exponential 
growth nature of the solution to flU). If we write the solution as 

X- = Acos(u k t), (6) 

where A is a slowly time dependent, complex function, the energy density can be 
given as follows: 

,2 !, ,2 

Pk = -\Xk\ + -\x-\ , 

= ^\A\ 2 u 2 k cos 2 (u k t) + - 1 A\W k sin 2 (uo k t) , 
= \\A?ul (7) 
From the equation @, we get 

Pk = ^2o-0 o m (/) cos(m (/) t)|y4| 2 cos 2 (u; fc t), 

= ^a<f> m^\A\ 2 cos(m^t)[l + cos(2u k t)], (8) 

and taking the average over one period of oscillations of inflaton field 0, the rhs of 
this equation will not vanish when the resonance condition uj k = holds so that, 

1 



Pk - -crfom^A 



2 



-Pk, (9) 



2 Analyzing in the narrow resonance region, we can regard that the amplitude of inflaton field 
is approximately constant. 
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where (|7|) is used in the second line. 

p k x can be identified as the production rate of the x-partide with momentum k 
due to parametric resonant decay of inflaton field <fi. Since the production rate is 
proportional to the energy density of produced x-P ar ticl es ( hence its occupation 
numbers in each k mode), we can see that this is an induced effect. That is, if 
many particles are created in the particular phase space, further production will be 
"induced" . 

Now we will find and use the explicit form of the solution for (|TJ) to obtain the 
previous result. Hamiltonian is a diagonalized operator for the free field, but in 
the presence of interaction, it is not diagonalized in terms of a and (V . It can be 
diagonalized at any instant of time by means of Bogolyubov transformations, which 
give the transition from the free field creation and annihilation operator, and a, to 
time dependent operators, b'(t) and b(t) (which preserve the commutation relation), 

b(t) = a{t)a + P*{t)a\ 

b\t) = P(t)a + a*(t)a\ (10) 
and a(t) and (3{t) can be written as [§ 

e « J VL k dt 
—i J Qkdt 

P = ^ (fi fc x- -**-), (ii) 

where fll = oj\ + 2a<f>, and the initial conditions for a, P are 

|a(0)| = l, /3(0) = 0. (12) 
Then the form of the solution can be written as, to the lowest order |J, 

X- = l— l x cos(^ fc t) — Ysm(u k t)), (13) 



where X, Y are complex weakly time-dependent functions. From Eqs. (|TT|) and (|12D 
the initial conditions for x become 

iX-(O) = w fc X-(0), |X-(0)| = (14) 
which lead to the initial values for X and Y: 

X(0)=e 4<5 , Y(0)=ie iS , (15) 
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where S is some phase which is not important and is set to zero below. 
We are only interested in the case where the resonance condition holds: 

co k - ^. (16) 

Therefore, the approximate solutions for X and Y are given by 

X ~ 1 + 1 4- 1 ~ - p-vt 
2 2' 

Y ~ i±le^*-— e-"*, (17) 



where /i = m^ 1 y o 2 <\>\ — A 2 , and A = uj"1 — Therefore we obtain the squares 

of the amplitude of X- and X- m t ne form 

\X-\ 2 = — [cosh(2^t) -sinh(2^t)sin(2a; fe t)], (18) 

|X-| 2 = — u^[cosh(2/4t) +sinh(2//t)sin(2a; fc t)], (19) 
2a> fc 

Taking the average over one period of oscillations of inflaton field and using 
Eqs. (|I~8| ) and ([19]), we obtain the averaged p£ and p£: 

(pj) ~ — cosh(2At) ~ ^ cosh(2At) (20) 
2 4 

(p£) ~ (T0o sinh(2pt)(sin(m^/f:) sin(2u; fc t)) = -<T0o sinh(2pt), (21) 

where the non-vanishing value can be achieved only when the resonance condition 
(16|) holds. Therefore, the following relation is found to be satisfied: 



p\ = — - tanh(2pt)pt ~ — -ft (as fit : large). (22) 

777.0 rricf, 



Here we get the same result as @. 

The total production rate is obtained by summing up all the mode: 

2tt 2 

' dt 



f uJi 1 f^ k 
J — du k -cr<f) sinh(2 J jj, 



duo' k 



dw' k ) 



2vr 2 I 2 J 2 C ° Sl1 m 3 



cosh rr I" , (23) 
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where p is regarded as adiabatically changing variable in the resonance band and 
the integration is taken over only inside the resonance band in the second line and, 
in the third, we use the resonance condition ( uo^ ~ ) and H is the Hubble's 

constant. Using the relation p = F^ ea 'p^, ( energy conservation ), we can calculate 
the resonant decay rate given by 



2n 2 p^ V 2 / 2 V Hml 



This is identical to the result of ref. (when J pdt is large). 

As mentioned before, we can obtain the similar results for other forms of inter- 
action between the inflaton and the x-particle or for |0 4 potential. [] Here, we only 
show the results corresponding to (|22|) . For the interaction £ int = —g4> 2 x 2 i we get 

« - < 25 > 

and for the theory V = j(p 4 , we obtain 

p\ ~ -^=pj, for £^ = -a0x 2 , (26) 



'° p* for £ m4 = -<70Y, (27) 



2c\/A 



where c ~ C(l). 

3 Restriction due to dissipative effects 

As shown in the previous section, the parametric resonant decay of inflaton field is an 
induced effect of produced bosonic particles. If there are large occupation numbers 
of produced bosons, further creation proportional to their occupation numbers will 
occur. The crucial point is that a lot of bosons can occupy in the same phase space. 

On the other hand, if particles in particular phase space decay into other particles 
or collide with each other, occupation numbers of bosons in that particular phase 
space cannot grow very fast. As a result, the resonant decay will be much suppressed. 



3 The behavior of inflaton field can be approximated with a sinusoidal function since its ampli- 
tude is small in the narrow resonance region. 
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In other words, the resonance band becomes narrower in the presence of decay or 
scattering processes. 

We begin with the mode equation for the x field which is produced by the decay 
of inflaton field 0. One convenient way to take into account the decay of the field \ 
is to introduce a dumping term in the mode equation. 4 Then the equation becomes 
a type of equation for damping oscillator with an external oscillating force for the 
theory V = ^m^cf) 2 and C int = —crcfix 2 (we can obtain the similar results for different 
types of effective potential V{4>) or interaction Lagrangian Ci nt ), 

Xk + F d Xk + [Wfc + 2cr0o sin(m 9i t)]xfc = 0, (28) 

where Td is the decay rate for the x field. Then we obtain the solution Xk for Eq.(p8|) 
in the following form: 

p ^ 

|X-| 2 = [cosh(2^) - sinh(2//*) sxa{2u k t)\. (29) 

The evolution of the energy density of the x-particle in k mode becomes 

d n 
di \2 



■ k 

Px 



ix-r+^ix-r) 



~^\X \ 2 -T d \x I 2 
2 dt lX - 1 dlX - 1 



- -p k x -T dP k x , (30) 

where we use Eq.([28|) in the second line and, in the third, take the average over 
one period of oscillations of inflaton field (p. The first term is the usual resonant 
production from decay of inflaton field <ft and the second is the energy decrease due 
to the decay of the x-particle. Therefore, the energy density of the x-particle in the 
mode k is reduced by the factor e _rd * compared with the case of no dissipation, and 
hence the considerable decrease of the resonant decay rate of inflaton field. 

Summing over mode k, we obtain the total resonant production of energy density 
of the x-particle: 

* = J 2^2^-1 



To include the dissipative process in the mode equation is not so simple as given in this 
paper. The precise treatment of dissipation is very complicated subject and beyond the scope of 
the present paper. However, we believe that the simple argument given here is reasonable in the 
physical point of view. 
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dwk-crcfro exp( 



27r 2 ~" re 4 



H 



(m±y 1 
2tt 2 V 2 / 4 



cxp 



/ 

V 



eft 



sm 



\ 



2(T</)n 



2a0 o \ 



1 - 



, 2(T0o 



(31) 



where, in the third line, we take into account that the resonance band width gets 

narrower due to dissipation: A = er0 o — > a(pQ\jl — ( ^g^ ) ^ ^ d < ' T( l ) o/ m <t>(— A 4 )' 
the production rate becomes approximately, 



f m<t> 
2vr 2 V 2 



exp 



27ra 2 2 
HrrA 



1 



2m0r d 



(32) 



Now we can see that the resonant production of x-particles is reduced due to dissi- 
pative effects, and hence the reduction of resonant decay rate of inflaton field <fi. 

Simpler explanation can be given as follows. The naive integration in the expo- 
nent in Eq.(|3~l~D becomes 



J (2/2 - Y d )dt = J 2fidt - J T d 
~ J 2fjdt - T d r re 



(33) 



where we regard T d as constant, and r res is the duration of the parametric resonance 
(the time while staying in the resonance band). If r res > T^ 1 , i.e., the lifetime of the 
X-particle is shorter than the resonant duration, the particle does not stay in the 
particular phase space long enough that the occupation number in that particular 
mode cannot grow so fast, which results in less resonant production. 



4 Conclusion 

We investigate the parametric resonant decay of inflaton field <fr and the dissipative 
effects of particles produced after the chaotic inflation. The equation for the particles 
produced can be described by Mathieu type equation, which has instability solution 
in some regions of parameters (instability bands). If the frequencies stay in the 
instability bands long enough, the solution will explosively grow, so that the number 
of created particles becomes exponentially large. Therefore the previous authors || 
H [lO], |ll| [L2| regarded that the parametric resonant decay is important. 
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In this article, we considered this phenomena from another point of view. The 
resonant decay is an induced effect of bosonic particles. If there are large occupation 
numbers of produced bosons, further creation proportional to their numbers will 
occur. A lot of bosons can occupy in the same phase space while one fermion can 
occupy in the particular phase space. This is the reason that the inflaton field does 
not decay explosively into fermions (Pauli's exclusion principle) [fTl[ P, p|, |T0[ . 



Then, in the case of decaying into bosons, if the dissipative effects prevent the oc- 
cupation numbers of bosons from growing very fast, the resonant decay will be much 
suppressed. The dissipative processes involve collisions between particles created, 
and their decay or annihilation into other lighter particles. 

The simplest way to look for this effect is compare the duration of resonance 
(the time that resonance frequency stays in the instability band) with the lifetime 
of produced particles. If the resonant duration is shorter than the particle lifetime, 
resonance occurs considerably. In the opposite case, parametric resonance should 
be extremely suppressed. We believe that this criterion is general. 

The dissipative effects are taken into account phenomenologically by TdXk term 
in the equation for the mode k of the field \i although it is difficult to obtain the 
term T d Xk from the first principle. Besides, decay rates of particles depend on its 
energy which is redshifted by the expansion of the Universe, so that the investigation 
is much more complicated. However, we expect that the result will remain the same 
qualitatively. 

There are some cosmo logical implication from this effects. First, those particles 
that are created by the decay of inflaton field should weakly interact with themselves 
or with other particles. Otherwise, the dissipative processes become effective so as 
to reduce the resonant decay of inflaton field considerably. Second, the reheating 
temperature will not be so high because thermalization takes for a long time. 

Finally we must mention in closing that our simple analytical investigation in this 
article for both parametric resonance and dissipative effects will be limited to only 
in the narrow resonance region where the effects of back-reaction to the amplitude of 
inflaton filed is rather small. There are some papers studying broad resonance, which 
is rather complicated f| [10, 11], 15]]. However, we expect that the basic physical 



process for broad resonance is the same as narrow one. Therefore, an induced 
nature of parametric resonant decay will also hold qualitatively in the broad resonace 
region so that the dissipation tends to suppress the resonant decay in the same way. 
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(It is model-dependent whether or not this suppresion becomes important, and we 
expect that it will be less efficient in the broad resonance region, since the resonance 
duration will be much shorter than the lifetime of created particles.) 

The authors are grateful to T. Yanagida for useful discussions. 
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